Evidence for new regimes in channeled granular flows 
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By means of extensive numerical simulations, we discovered new regimes in channeled granular 
flows. These regimes emerge from the destabilization of unidirectional flows upon increase of mass 
holdup and inclination angle. They are characterized by no n- trivial internal structures, with sec- 
ondary flows and/or an heterogeneous mass distribution. The newly identified states include flows 
with a pair of rolls, superposed rolls, a basal ordered layer topped by secondary flows, and states 
where a dense core is supported by a gaseous phase. Besides, the supported regimes exhibit striking 
behaviors at high mass holdup such as symmetry breaking leading to periodic oscillations of the 
dense core. Interestingly, these different regimes are shown to obey a universal scaling law for the 
mass flow rate as a function of the mass hold-up. These new states provide a unique set of 3D flow 
regimes on which granular rheological models can be tested and open perspectives for interpreting 
geophysical granular flows. 

PACS numbers: 47.57.Gc, 45.70.-n 



Introduction- Granular flows are encountered in 
many geophysical flows and industrial processes but there 
is still no general theoretical description for these flows. 
Progress has been made for flows down inclines [T]-[3] and 
current theories give a reasonable description for unidi- 
rectional and uniform flowing states. These simple flows 
have been shown to be unstable at large inclination an- 
gles and to develop longitudinal vortices [11 15]. However, 
little is known about the instabilities of these unidirec- 
tional flows. One objective reason for this is that the pa- 
rameter space accessible to experiments and simulations 
is restricted to a reduced domain compared to what is 
possible in nature. For example, the experimental obser- 
vations of steady and fully developed flows at high an- 
gle and high mass flow rate would require long channels. 
We have developed optimized discrete element method 
(DEM) simulations that allow exploration of granular 
flows in wide channels at high inclination angles. 

We consider a simple flow geometry with flat frictional 
basal and lateral walls. This geometry has been already 
used in previous experimental [6] and numerical [7] works 
and includes a wide variety of industrial or geophysical 
flows. Our numerical simulations reveal the existence of 
a rich variety of new regimes never reported in the litera- 
ture. These new regimes emerge from the destabilization 
of steady and fully developed unidirectional flows upon 
increase of the mass holdup and the slope. They exhibit 
non-trivial internal structure with secondary flows and 
heterogeneous mass distribution, including states with 
array of rolls, a pair of rolls atop a basal sheared ordered 
layer, and dense cores supported by a gaseous phase. 
These new flow regimes are expected to be encountered 
in many real situations and raise a new challenge for ex- 
tending the scope of current rheological models [2j. 

Method- We have performed 3D simulations of gran- 
ular flows with a classical discrete element model fT.'S^.'Qj. 



This technique consists in integrating the equations of 
motion of all grains, and resolving their interactions by 
considering them as soft interpenetrating spheres. De- 
tails about this method can be found in our previous 
work [7J . The material parameters were chosen to match 
the experiments of Louge and Keast grain density 
pg = 2.5 X lO^kg/m^, grain diameter D = 3mm^ coef- 
ficient of restitution e = 0.972, grain/grain friction co- 
efficient fig = 0.33 and grain/wall friction fi^n = 0.59. 
We used the same channel geometry as that employed 
by Louge and Keast in their experiments [6j. The chan- 
nel has a rectangular section and a width W ^ 68D. 
We explicitly account for the presence of lateral walls, 
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Figure 1: Granular flows down inclined channels phase dia- 
gram. Circled letters are used throughout this document to 
identify the regimes. Each of the 385 gray dots represents the 
result of a simulation that has reached its final state. 
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Figure 2: (Color online). Snapshots of different regimes. Packing fraction (a) and velocity field (b) in the transverse plane 
for the rolls regime @. Packing fraction for the transition (c) to the supported symmetric dense core (d) regime (c). Packing 
fraction for the supported regime at high mass holdup ®, presenting an oscillating asymmetry (e). (f) Average velocity field 
in the transverse plane showing stacked rolls @ at high mass holdups. Average packing fraction (g) and velocity (h) in the 
transverse plane, presenting rolls in a "melted" zone on top of a base consisting of sheared ordered layers (b). 



unlike most numerical studies considering periodic con- 
ditions in the flow transverse direction. Walls are not 
only physically relevant, but necessary in order to reach 
steady flows at high angles 6 in simulations. They are 
sufficiently far away to allow for the development of sec- 
ondary flows, which would not appear in more confined 
setups. Simulations were started with densely packed 
grains dropped at a small altitude from the basal wall 
and were run up to reaching a steady regime, possibly 
oscillatory. 

We have varied the angle of inclination and the 
amount of matter per unit surface measured via the mass 
holdup H defined as: H = ^ u{z)dz^ where u is the 
particle volume fraction. H is the dimensionless height 
that the flow would have without voids between grains. 

Phase Diagram- Qualitatively distinct regimes are 
identified according to the values of the inclination angle 
and mass holdup (see Fig. [T] and |2|. Larger figures and 
snapshots are available in the supplementary materials. 
Fig. [l] indicates in the parameter space (^, H) the domain 
of existence of the different identified regimes. These are 
labeled by circled letters and are briefly described below: 

- Regime (u) corresponds to Unidirectional flows with 
sheared ordered layers of grains. The properties of these 
flows match the experimental observations reported by 
Louge and Keast [6j as detailed in [7J. 

- Regime @ corresponds to the presence of secondary 
Rolls, visible in Fig. These structures were previ- 
ously reported and analyzed in recent experimental and 
numerical works [4l|5][7]. In particular, an analogy with 
Rayleigh-Benard convection rolls in heated liquid layers 
can be made, where the "granular temperature" plays the 
same role as the classical thermodynamics temperature 
[4J. Rolls break at larger angles when the flow lifts up 
on top of a layer of granular gas (regime (c)), as illus- 
trated in Fig. 



- Regime (c) is characterized by the presence of a dense 
Core floating on top of a granular gaseous phase. It is 
clearly visible on the density map of Fig. [2]i. This regime 
is referred to as a supported flow. The core slowly "evap- 
orates" as the angle increases for a gradual transition to 
granular gas at larger angles (shown in the supplemen- 
tary materials). When the mass holdup increases, the 
core starts to swing back and forth from left to right and 
loses its axial symmetry. 

- Regime (a) denotes the Asymmetric core. For larger 
H and 6>, a "plume" forms on top of the core, as shown 
in Fig. [2^. This structure oscillates from one side to the 
other through time, recovering the symmetry on average. 

- Regime (s) corresponds to the Superposed rolls and 
appears at larger mass holdups H than regime (r). An 
example is shown in Fig. [2]F. 

- Regime (b) is characterized by the presence of a basal 
layered structure. The observed order (see Fig. [2^) is dy- 
namically maintained by collisions and cage effects. The 
layers are sheared and not static. Rolls are present in the 
disordered zone on the top of the basal layers and are lo- 
calized close to the lateral walls (Fig. [2]i). This complex 
flow regime that exhibits both order /disorder transition 
and secondary flows raises a significant challenge for the- 
oretical modeling. 

These different flow regimes open many perspectives 
given the generality of the inclined channel geometry. 
The regimes (u) and (E) have been analyzed extensively in 
our preceding work [7J. Here, we focus on the supported 
regimes and present further details below. 

Supported regimes- These regimes are of great inter- 
est because they may explain part of why low friction 
is observed for long run-out landslides. Campbell [lOj 
indeed suggested that the existence of a layer of highly 
agitated particles at low concentration beneath a densely 
packed main body could reduce the apparent basal fric- 
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tion. By employing 2D [TT] and 3D [12] discrete particles 
numerical simulations, it has been shown that a granu- 
lar material subjected to gravity approaches this state if 
its initial kinetic energy is large enough. All these flows 
are in general not steady but decelerate or accelerate de- 
pending on the slope. Taberlet et al [12] simulated flows 
on bumpy inclined with lateral periodic boundary con- 
ditions, but these supported flows reach a steady regime 
only for a unique value of the slope. In our case, steady 
supported flows are obtained in a large range of incli- 
nation angles. This regime is very stable and robust to 
variations of the model parameters. 

Fig. [3] shows typical density profiles for the supported 
regimes. Regimes (a) and @ both present a dense core 
floating above a granular gaseous phase (Fig. [2]l,e), and 
toped by a dilute "atmosphere". It is possible to see the 
signature of the "plume" (regime A) in the profiles above 
the core. The center of mass of the flow is located just 
on top of the core, so a large fraction of the matter is 
spread in the dilute atmosphere over large distance. 

When the mass holdup increases the core lifts up and 
densifies, as shown in Fig. [3^. The lateral width of the 
core (see Fig. |2]i,e) decreases with H because the lateral 
pressure pushes the grains toward the central core. This 
core can reach very high values of the volume fraction 
up to 0.6 at large mass holdup, while the density in the 
supporting basal gaseous layer is below 0.2. The volume 
fraction in the dilute atmosphere above the core is well 
fitted by a decreasing exponential: u{z) ex exp(— z/i^c), 
where He represents the characteristic height of the at- 
mosphere (see the supplementary materials). 

Surprisingly, when the angle increases, the average 
altitude of the core remains constant (see Fig.[3]3). How- 
ever, the center of mass lifts up and the core thickness 
decreases as a non-negligible part of the material is trans- 
ferred into the top granular gaseous phase. Fig. shows 
that for the supported regimes (indicated by dots in the 
Figure), the position Cm of the center of mass increases 
linearly with tan^: Cm = atan^ + 6, where the slope 
a is independent of the mass holdup. In contrast, the 
parameter b increases with increasing mass holdup and 
reflects the correspondingincrease of the core thickness 
with H described in Fig. [3p. 

The vertical expansion of the flow is necessary to in- 
crease its friction on the lateral boundaries and to balance 
the driving force which increases with the inclination an- 
gle. The basal friction cannot exceed /j^mMg cos 6, where 
/irn = 0.59 is the microscopic friction value used in the 
simulations [7J and M is the mass of the grains. Thus 
for large angles, a large part of the friction comes from 
the lateral walls. The apparent friction coefficients on 
the base /i5 and walls /i^^ are computed as the ratio of 
tangential to normal stresses averaged over all simulated 
impacts. Fig. shows the dependency of /i5 on and 
H. ji^ follows the same trend up to very minor varia- 
tions (see inset of Fig.|4]3). Both coefficients increase and 




0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 



Average density, channel center ±5D 

Figure 3: (Color online). Vertical profiles of the density of the 
flow at the center of the channel, averaged over lOD. (a) for 
a fixed angle = 42° when varying the mass holdup, and (b) 
for a fixed mass holdup ^=8 when varying the angle. Empty 
dots and triangles indicate the position of the center of mass. 

saturate at high inclination angles. In contrast, they de- 
crease with the mass holdup: for a given angle, the basal 
friction reduces as more matter is added in the flow. This 
effect is reminiscent of what is observed for large natural 
flows in long run-out avalanches [lOj. 

Using a force balance on the periodic cell in the station- 
ary regime, we write: Mg sin 6 = 2/i^L Nyj{z)dz + 
/JjIjNijWL, with Nw{z) and the normal stress at side- 
walls and at the base respectively. L = 20D is the 
length of the cell in the flow direction. The first term 
is the friction of the wall on the flow, and the second 
the friction of the base. We assumed Coulomb's friction 
as both /j^w and /i^ saturate in the supported regimes 
(Fig. ^p)' The force balance in the other direction gives: 
MgcosO = NijWL. This simple model thus predicts: 
2 N^{z)dz/Nb = ^ (tan (9 - //?,). The first term rep- 
resents a characteristic height of the flow and will be 
denoted by Hp (if the pressure were purely hydrostatic. 
Hp would correspond exactly to the height of the flow). 
This relation expresses the fact that the characteristic 
height of the flow Hp varies linearly with tan^. This ex- 
plains why the position Cm of the center of mass, which 
is of the order of the half height of the flow, evolves 
linearly with tan^. Remarkably, we find here a slope 
a = 54.6 ~ W/ (2/i^) in agreement with the hypothesis 
that Cm ~ Hp/2. We also find that the "atmospheric 
height" He increases linearly with tan in the supported 
regimes. The slope here is about two third of Wf/j.^, 
corresponding to He ~ 3Hp/2. 

Transient- The typical temporal evolution of the av- 
erage flow velocity to reach the steady regime is shown 
in Fig. |5|i for ^ = 5. The data can be fitted with 
good accuracy (RMS error less than 6%) by an expo- 
nential saturation toward the limit velocity: V{t) = 
Vl — {Vl — Vq) exjp {—t / r) ^ where Vq is the initial aver- 
age velocity in the flow direction, Vl the limit velocity 
and r a characteristic time. The adaptation of the flow 
during the transient requires a strong reorganization of 
its structure. Outstandingly, in spite of this, we observe 
a very regular evolution, with a similar behavior (expo- 
nential saturation) for all the angles, and thus for all the 
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Figure 4: (Color online), (a) The center of mass Cm altitude, 
with respect to the slope angle 0. Markers indicate the config- 
urations in the supported regimes @ and (c) and (b) effective 
friction coefficients on the base fih and the wall {fiw^ inset), 
as defined in the main text. 
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Figure 5: (Color online), (a) Typical temporal evolution of 
the mean flow velocity, for i7 = 5 and ^ = 12° to ^ = 48° 
every 3°. (b) Characteristic time for the flow to reach the limit 
velocity (small angles and intermittent states not shown). 



regimes. On the other hand, the characteristic time r 
is affected by changes of the inclination angle 6 and the 
mass holdup as can be seen in Fig.jsjD, which reveals a 
general tendency for longer transients in the roll regimes 
(r) and (s) in comparison with supported regimes (c) and 

(a) at larger angles. 

The exponential velocity saturation suggests that the 
flow experiences a viscous-like drag force proportional to 
the velocity: Fy = —MV/r. Such a viscous frictional 
force due to collisions of grains with boundaries was al- 
ready reported a long time ago [13j. The viscous damp- 
ing constant a = M/r follows the inverse trend of r, 
and increases with the mass holdup (shown in the sup- 
plementary materials), a is large for the regimes where 
grain interactions are dominated by collisions: (u) and 

(b) where collisions maintain dynamically the order, and 
supported regimes (c) and (a) lubricated by a highly col- 
lisional basal gaseous layer. 

Steady state velocity- In contrast with the complex 
dependency of the relaxation time r on ^ and we 
identify a simple universal dependency for the limit ve- 
locity: Vl/H'^^^ ^ A^iiiO + 5, where A and B are con- 
stant coefficients. In Fig. [6) we plotted the limit velocity 
Vl rescaled by H'^^^ as a function of the inclination an- 
gle for various mass holdup. A remarkable collapse is 
observed, given the large diversity of the flow regimes. 
In steady state, the mass flow rate is simply given by 
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Figure 6: (Color online). Scaling law for the steady state 
velocity: Vl/H^^^ ^ AsinO + B, with A ^ 122 and B ^ -37. 
The collapse is not perfect, but nonetheless remarkable given 
the diversity of regimes. 

Q = VlH. Our results thus show that Q cx H^^^ for a 
fixed angle ^, which is close to the scaling law (Q oc H^^'^) 
derived by Louge and Keast from experiments [6j but ob- 
tained over a smaller range of flow rate Q. 

Conclusion- We explore granular flows configura- 
tions that were not previously attained in experimental 
or numerical works, with wide channels and large incli- 
nation angles. Using a simple flow configuration with 
flat lateral and basal boundaries, we have discovered, by 
varying the inclination angle and mass holdup, a mul- 
titude of new regimes ranging from rolls to dense cores 
supported by a highly agitated granular gaseous phase. 
These regimes present non-trivial features including het- 
erogenous density, secondary flows, anisotropy, symme- 
try breaking, and dynamically maintained order. 

For increasing angles, the flow inflates to balance the 
growing driving force by a larger friction on the lateral 
walls. At some stage the bulk of material lifts up to ride 
on a thin layer of agitated grains (supported regimes). 
For increasing mass holdup the lateral and basal ef- 
fective friction coefficients reduce, revealing properties 
similar to long run-out avalanches in nature. Besides, we 
showed that the mass flow rate of these flows obeys a uni- 
versal scaling law in terms of H. The new identified states 
provide a unique set of granular flow regimes for testing 
3D rheological models and making further progress in 
the modeling of granular flows. These results also open 
new perspectives for understanding geophysical flows and 
should encourage investigations of granular flows in wider 
channels and for higher mass holdup. 
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